Particle Production During Inflation: Observational Constraints and Signatures 
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In a variety of inflation models the motion of the inflaton may trigger the production of some 
non-inflaton particles during inflation, for example via parametric resonance or a phase transition. 
Particle production during inflation leads to observables in the cosmological fluctuations, such as 
features in the primordial power spectrum and also nongaussianities. Here we focus on a prototype 
scenario with inflaton, 4>, and iso-inflaton, y, fields interacting during inflation via the coupling 
g 2 (4>~ 0o) 2 X 2 - Since several previous investigations have hinted at the presence of localized "glitches" 
in the observed primordial power spectrum, which are inconsistent with the simplest power-law 
model, it is interesting to determine the extent to which such anomalies can be explained by this 
simple and microscopically well-motivated inflation model. Our prototype scenario predicts a bump- 
like feature in the primordial power spectrum, rather than an oscillatory "ringing" pattern as has 
previously been assumed. We discuss the observational constraints on such features using a variety 
of cosmological data sets. We find that bumps with amplitude as large as 0(10%) of the usual 
scale invariant fluctuations from inflation, corresponding to g 2 ~ 0.01, are allowed on scales relevant 
for Cosmic Microwave Background experiments. Our results imply an upper limit on the coupling 
g 2 (for a given 4>o) which is crucial for assessing the detectability of the nongaussianity produced 
by inflationary particle production. We also discuss more complicated features that result from 
superposing multiple instances of particle production. Finally, we point to a number of microscopic 
realizations of this scenario in string theory and supersymmetry and discuss the implications of our 
constraints for the popular brane/axion monodromy inflation models. 

PACS numbers: 11.25.Wx, 98.80.Cq 



I. INTRODUCTION 

Recently, there has been considerable interest in in- 
flationary models where the motion of the inflaton trig- 
gers the production of some non- inflation (iso-curvature) 
particles during inflation 0, 0, i, i, [S IS S B ©, EE 
El EE 0, 0, EE EE E3- Examples have been stud- 
ied where this particle production occurs via parametric 
resonance |l|, [2|, [3|. |j, l a , |fj [ , as a result of a phase tran- 
sition 0, HTI1, EE TUT EE EH, or otherwise. In some 
scenarios, backreaction effects from particle production 
can slow the motion of the inflaton on a steep poten- 
tial [3 EE EH) providing a new inflationary mechanism. 
Moreover, inflationary particle production arises natu- 
rally in a number of realistic microscopic models from 
string theory [MEE EE EE El, HI and also supersym- 
metry (SUSY) [21(~^ 

In [5| it was shown that the production of massive iso- 
curvature particles during inflation (and their subsequent 
interactions with the slow roll condensate) provides a 
qualitatively new mechanism for generating cosmological 
perturbations. This new mechanism leads to a variety of 
novel observable signatures, such as features [j| and non- 
gaussianities [Ellli m the primordial fluctuations. In this 
paper we study in detail the observational constraints on 
such distortions of the primordial power spectrum for a 
variety of scenarios. 

One motivation for this study is to determine whether 
features generated by particle production during in- 
flation can explain some of the anomalies in the ob- 



served primordial power spectrum, P(k). A number 
of different studies have hinted at the possible pres- 
ence of some localized features in the power s pec trum 

[I El, HE Hi HE HE Hi HE Hi, HE M HE HE 11, 

which are not compatible with the simplest power law 
P(k) ~ fc™ a_1 model. Although such glitches may sim- 
ply be statistical anomalies, there is also the tantalizing 
possibility that they represent a signature of primordial 
physics beyond the simplest slow roll inflation scenario. 
Forthcoming polarization data may play a crucial role in 
distinguishing between these possibilities 13]. However, 
in the meantime, it is interesting to determine the extent 
to which such features may be explained by microscopi- 
cally realistic inflation models. 

We consider a very simple model where the inflaton, 
0, and iso-inflaton, x, fields interact via the coupling 
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(1) 



We focus on this simple prototype model in order to il- 
lustrate the basic phenomenology of particle production 
during inflation, however, we expect our results to gener- 
alize in a straightforward way to more complicated sce- 
narios. Models of the type ([1]) have been considered as a 
probe of Planck-scale effects [l| and offer a novel exam- 
ple of the non-decoupling of high energy physics during 
inflation. 1 



For reasonable values of g the % particles are extremely mas- 
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At the moment when cf> = <^>o (which we assume occurs 
during the observable range of e-foldings of inflation) the 
X particles become instantaneously massless and are pro- 
duced by quantum effects. This burst of particle produc- 
tion drains energy from the condensate 4>(t), temporarily 
slowing the motion of the inflaton background and vio- 
lating slow roll. Shortly after this moment the x particles 
become extremely non-relativistic, so that their number 
density dilutes as a -3 , and eventually the inflaton settles 
back onto the slow roll attractor. 

Several previous papers [2, H, H, H[ have studied the 
temporary slowing-down of the inflaton background us- 
ing the mean-field equation 

4> + 3Hj> + V^+g 2 (0-fo)(x 2 ) =0 (2) 

where the vacuum average (x 2 ) is computed following 
[3rj l37l ]. Using this approach one finds that the tran- 
sient violation of slow roll leads to a "ringing pattern" 
(damped oscillations) in the spectrum of cosmological 
fluctuations leaving the horizon near the moment when 
4> = 4>o 01 • I n @j U observational constraints on particle 
production during inflation were discussed in the context 
of this mean field treatment. 

However, in [5] cosmological fluctuations in the model 
(fTj) were re-considered, going beyond the mean-field 
treatment of <p. 2 It was pointed out that the massive 
X particles can rescatter off the condensate to generate 
bremsstrahlung radiation of long-wavelength 8<p fluctua- 
tions via diagrams such as Fig. [TJ Multiple rescattering 
processes lead to a rapid cascade of power into the infra- 
red (IR) - IR cascading. The inflaton modes generated 
by IR cascading freeze once their wavelength crosses the 
horizon and lead to a bump-like feature in the cosmologi- 
cal perturbations that is illustrated in Fig. [5] This feature 
is complimentary to the usual (nearly) scale-invariant 
quantum vacuum fluctuations from inflation. The bump 
dominates over the ringing pattern discussed above by 
many orders of magnitude, independently of the value of 
9 2 - 



x 




FIG. 1: Rescattering diagram. 



sive for almost the entire duration of inflation excepting a tiny 
interval, much less than an e-folding, about the point <f> = <f>0- 
However, the \ field cannot be integrated out due to the non- 
adiabatic time dependence of the mode functions, see |35|| for 
further discussion. 
2 See also [l5l | for a complimentary analysis and 38] for a detailed 
analytical treatment of the dynamics. 
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FIG. 2: The bump-like features generated by IR cascading. 
We plot the feature power spectrum obtained from fully non- 
linear lattice field theory simulations (the red points) and also 
the result of an analytical calculation (the dashed blue curve) 
using the formalism described in |38j]. We also superpose the 
fitting function ~ k 3 e~ 7vk /(2fe * ) (the solid black curve) to il- 
lustrate the accuracy of this simple formula 

In light of the results of [B[ it is clear that the obser- 
vational constraints on the model |T]) need to be recon- 
sidered. Since previous studies have suggested marginal 
evidence for localized power excesses in the CMB using 
both parametric @, and non-parametric [3(| HH tech- 
niques, it is interesting to determine if a simple and well- 
motivated model such as ([!} can explain these anomalies. 
To answer this question we provide a simple semi-analytic 
fitting function that accurately captures the shape of 
the feature generated by particle production and IR cas- 
cading during inflation. Next, we confront this modi- 
fied power spectrum with a variety of observational data 
sets. We find no evidence for a detection, however, we 
note that observations are consistent with relatively large 
spectral distortions of the type predicted by the model 
([1]) . If the feature is located on scales relevant for Cos- 
mic Microwave Background (CMB) experiments then its 
amplitude may be as large as 0(10%) of the usual scale- 
invariant fluctuations, corresponding to g 2 ~ 0.01. Our 
results translate into a ^-dependent bound on the cou- 
pling g 2 which is crucial in order to determine whether 
the nongaussian signal associated with particle produc- 
tion and IR cascading is detectable in future missions 

We also consider the more complicated features which 
result from multiple bursts of particle production and IR 
cascading. Such features are a prediction of a number 
of string theory inflation models, including brane/axion 
monodromy [HI, [lfj [2(| ■ For appropriate choice of the 
spacing between the features, we find that the constraint 
on g 2 in this scenario is even weaker than the single-bump 
case. 

Although we focus on the interaction |T]), our findings 
may have some bearing also on models with phase transi- 
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tions during inflation. A simple prototype model for the 
latter scenario is 

v{fa x ) = \{x 2 -v 2 ) 2 + 9 ^x 2 tf + v i M 

At the moment when cf> — y/\v / g massive iso-inflaton 
X particles are produced copiously by tachyonic (spin- 
odal) instability (39|. These produced particles will sub- 
sequently interact with the condensate via rescattering 
diagrams similar to Fig. [I] Hence, we expect the features 
produced by inflationary phase transitions to be quali- 
tatively similar to the bumps considered in this paper. 
(This intuition is consistent with a second order com- 
putation of the cosmological perturbations in a closely 
related model • See (3J| for a discussion of the obser- 
vational consequences.) 

In the literature it is sometimes argued that inflation- 
ary phase transitions can be studied using a toy model 
with a sharp step-like feature in the inflaton potential. 
This potential-step model predicts a ringing pattern in 
the power spectrum, very much analogous to the mean 
field treatment of resonant particle production during in- 
flation, discussed above. This treatment does not take 
into account the violent growth of inhomogeneities of the 
fields that occurs during realistic phase transitions [39| 
and, in particular, does not capture rescattering effects 
and IR cascading. In the case of resonant particle pro- 
duction, these nonlinear effects have a huge impact on 
the cosmological fluctuations Hence, it is far from 
clear if the potential-step model provides a good effec- 
tive description of inflationary phase transitions. 3 

Of course, inflation models with steps in V(fa) (or 
its derivatives) may be considered on phenomenological 
grounds, irrespective of the motivation from inflationary 
phase transitions. In [4(| 5l[ cosmological perturbations 
from models with step-like features and discontinuities 
in higher derivatives were considered, as were the micro- 
scopic motivations for such constructions. See [13, |43[ 
for a study of the nongaussianities induced in a variety 
of single-field models with steps or oscillations in the in- 
flaton potential. 

The outline of this paper is as follows. In section |TT] 
we provide a simple parametrization of the features that 
are imprinted on the primordial power spectrum by one 
or more bursts of particle production during inflation. In 
section Hill we describe our method and discuss the obser- 
vational data sets employed to derive constraints on this 
modified power spectrum. In section IIVI we present ob- 
servational constraints on various scenarios. In section [V] 
we present some microscopic realizations of our scenario 
and discuss the implications of our findings for popular 
string theory/SUSY inflation models with a special em- 
phasis on brane monodromy. Finally, in section I VII we 
conclude. 



3 See also [44ll for a related discussion. 



II. A SIMPLE PARAMETRIZATION OF THE 
POWER SPECTRUM 

In Q it was shown that particle production and IR 
cascading during inflation in the model (TT]) generates 
a bump-like contribution to the primordial power spec- 
trum. As shown in Fig. [2J this feature can be fit with 
a very simple function fl, ump ~ k z e~^ k /( 2k *\ The 
bump-like contribution from IR cascading is complimen- 
tary to the usual (nearly) scale-invariant contribution to 
the primordial power spectrum P vac ~ fc" 3-1 coming 
from the quantum vacuum fluctuations of the inflaton. 
The total, observable, power spectrum in the model (TT]) 
is simply the superposition of these two contributions: 
P(k) ~ fc" s_1 + k 3 e~^ k /( 2fc *). This simple formula can 
be motivated from analytical considerations [D, |38| and 
provides a good fit to lattice field theory simulations near 
the peak of the feature and also in the IR tail. 4 

It is straightforward to generalize this discussion to 
allow for multiple bursts of particle production during 
inflation. Suppose there are multiple points cj> = 4>i 
(i = 1, ••• , n) along the inflationary trajectory where 
new degrees of freedom Xi become massless: 

71 n 2 

^t-EyW^^ (3) 

i=0 

For each instant tj when 4> = fa there will be an associ- 
ated burst of particle production and subsequent rescat- 
tering of the produced massive Xi on " the condensate fat). 
Each of these events proceeds as described above and 
leads to a new bump- like contribution to the power spec- 
trum. These features simply superpose owing to that fact 
that each field \i is statistically independent (so that the 
cross terms involving XiXj with i ^ j in the computation 
of the two-point function must vanish). Thus, we arrive 
at the following parametrization of the primordial power 
spectrum in models with particle production during in- 
flation: 




where A s is the amplitude of the usual nearly scale in- 
variant vacuum fluctuations from inflation and fco is the 
pivot, which we choose to be fco = 0.002 Mpc -1 following 
[4a ]. The constants Ai depend on the couplings gf and 
measure the size of the features from particle production. 
We have normalized these amplitudes so that the power 



4 This fitting formula does not capture the small oscillatory struc- 
ture in the UV tail of the feature (see Fig. |2j which does not 
concern us since that region is not phenomenologically interest- 
ing. 
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in the i-th bump, measured at the peak of the feature, is 
given by A- The location of each feature, k it is related 
to the number of e-foldings N from the end of inflation 
to the time when the i-th burst of particle production 
occurs: roughly ln(fcj/iJ) ~ Ni where N = JV» at the 
moment when <f> = fa. From a purely phenomenological 
perspective the locations ki are completely arbitrary. 

We compare ^ to lattice field theory simulations in 
order to determine the amplitude Aj in terms of model 
parameters. We find 

A-l.0ixicr 6 5l 15/4 (5) 

Assuming standard chaotic inflation V = m 2 fa 2 /2 we 
have tested this formula for g 2 — 1, 0.1, 0.01, taking both 
fa = 2M P and fa = 3.2M p . We found agreement up to 
factors order unity in all cases. 

Theoretical consistency of our calculation of the shape 
of the feature bounds the coupling as 10~ 7 ;$ gf ,$ 1 [j|. 
Hence, the power spectrum (jj) can be obtained from sen- 
sible microphysics only when 10~ 20 < A.- L <, 10 -6 . This 
constraint still allows for a huge range of observational 
possibilities: near the upper bound the feature is con- 
siderably larger than the vacuum fluctuations while near 
the lower bound the feature is completely undetectable. 

Note that for each bump in Q the IR tail Pbump — > k 3 
as k — > is similar to the feature considered by Hoi, 
Cline & Holder in [33], consistent with causality argu- 
ments about the generation of curvature perturbations 
by local physics. 



III. DATA SETS AND ANALYSIS 



We employ the Union Supernova la data (307 SN la 
samples) from The Supernova Cosmology Project 

Large Scale Structure (LSS) 

The 2dF Galaxy Redshift Survey (2dFGRS) data 
and Sloan Digital Sky Survey (SDSS) Lumi- 
nous Red Galaxy (LRG) data release 4 [68| are uti- 
lized. 

Note that we have used the likelihood code based 
on the non-linear modeling by Tegmark et al. [68j 
(marginalizing the bias b and the Q parameter). 
However with a large bump in the linear power 
spectrum, this naive treatment may be not suffi- 
cient to characterize the non-linear response to the 
feature on small scales. Ideally, this should be ob- 
tained from N-body simulations, however, such a 
study is beyond the scope of the current work. 

There are several other caveats on our results in 
the high-fc regime. First, we assume linear bias for 
the galaxies, which may not be entirely safe at suffi- 
ciently small scales. Moreover, sharp features in the 
matter power spectrum can cause sharp features in 
the bias as a function of k. 

Keeping in mind these caveats our constraints on 
small scales k £ 0.1 Mpc -1 should be taken with a 
grain of salt and considered as accurate only up to 
factors order unity. 

Weak Lensing (WL) 

Five WL data sets are used in this paper. The ef- 
fective survey area Aff and galaxy number density 
n e ff of each survey are listed in Table |TJ 



The primordial power spectrum for our model is 
parametrized as Q. Our aim is to derive observational 
constraints on the various model parameters A s , n s , ki 
and Aj using CMB, galaxy power spectrum and weak 
lensing data. To this end we use the cosmoMC package 
[11] to run Markov Chain Monte Carlo (MCMC) cal- 
culations to determine the likelihood of the cosmological 
parameters, including our new parameters A and fcj. We 
employ the following data sets. 



TABLE I: Weak Lensing Data Sets 



Data sets 


A e g 


n e ff 




(deg 2 ) 


(arcmin -2 ) 


COSMOS [69, 75] 


1.6 


40 


CFHTLS-wide [70, 71] 


22 


12 


GaBODS [72, 73] 


13 


12.5 


RCS [72, 73] 


53 


8 


VIRMOS-DESCART [71, 74] 


8.5 


15 



Cosmic Microwave Background (CMB) 

Ou r complete CMB data sets include WMAP-5yr 
[.BOOMERANG MJM, HI, ACBAR 
Hi, H|, CBI [H[5EjM,[5|, vsa [H, dasi 




|61J, and MAXIMA [6J|. We have included 
the Sunyaev-Zeldovic (SZ) secondary anisotropy 
[H [H for WMAP-5yr, ACBAR and CBI data 
sets. The SZ template is obtained from hydrody- 
namical simulation [65]. Also included for theoret- 
ical calculation of CMB power spectra is the CMB 
lensing contribution. 

Type la Supernova (SN) 



For COSMOS data we use the CosmoMC plug- 
in written by Julien Lesgourgues [75| . modified to 
do numerical marginalization on three nuisance pa- 
rameters in the original code. 

For the other four weak lensing data sets we use 
the likelihood given by [76^ . To calculate the like- 
lihood we have written a CosmoMC plug-in code, 
with simplified marginalization on the parameters 
of galaxy number density function n(z). More de- 
tails about this plug-in can be found in flT\ - 

As for the LSS data, for small scales k > 0.1 Mpc -1 
there is the caveat that the nonlinear evolution of 
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TABLE II: Constraints on the standard ("vanilla") cosmo- 
logical parameters for the single-bump model. All errors are 
95.4% confidence level 





0.0227t7 00 l 




n 1199 +U.UUbU 
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1.0424±™ 
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0.081^ 


n s 


u.yoo_ 024 


ln[10 1( M s ] 


3 206 +U UV4 


A S z 


1 62 +u - Y(i 




0.264ir ^ 






z re 


10.5+2:7 


H 


71.6+2:4 



the power spectrum in the presence of bump-like 
distortions may not be treated accurately. 

IV. OBSERVATIONAL CONSTRAINTS 

We now present our results for the observational con- 
straints on particle production during inflation, assuming 
two different scenarios. 



A. A Single Burst of Particle Production 

The minimal scenario to consider is a single burst of 
particle production during inflation, which corresponds 
to taking n = 1 in ([3]). The power spectrum is given 
by (j4]) with n = 1 and, with some abuse of notation, 
we denote k\ = km and A± = Am- The prior we have 
used for Am. is < Am < 25 x 1CT 10 , and for fci R is 
— 9.5 < ln[fc/Mpc _1 ] < 1. The former condition ensures 
that the bump-like feature from IR cascading does not 
dominate over the observed scale invariant fluctuations 
while the latter is necessary in order to have the feature 
in the observable range of scales. In Fig. [3] we plot the 
marginalized posterior likelihood for the new parameters 
Am and fciR describing the magnitude and location of 
the bump while in Table |TT] we give the best fit values for 
the remaining (vanilla) cosmological parameters. 

For very large scales <, Gpc -1 , the data do not contain 
much information (due to cosmic variance) and hence 
the constraint on any modification of the power spec- 
trum is weak. In this region the spectral distortion may 
be larger than 100% of the usual scale invariant fluc- 
tuations and couplings g 2 order unity are allowed. For 
smaller scales k > Gpc -1 the constraints are stronger 
and we have, very roughly, Am/ A s < 0.1 corresponding 
to g 2 < 0.01. For very small scales, k >, 0.1 Mpc -1 our 
constraints should be taken with a grain of salt since the 
nonlinear evolution of the power spectrum may not be 
modeled correctly in the presence of bump-like distor- 
tions. At small scales nonlinear effects tend to wipe out 



Nonlinear scales 



x 0.5 









Wk 


■ 


CMB WL + 


LSS 



FIG. 3: Marginalized posterior likelihood contours for the 
parameters Air and km (the magnitude and position of the 
feature, respectively) in the single-bump model. Black and 
grey regions correspond to parameter values allowed at 95.4% 
and 99.7% confidence levels, respectively. At small scales, to 
the right of the dashed vertical line, our results should be 
taken with a grain of salt since the nonlinear evolution of the 
power spectrum may not be modeled correctly in the presence 
of bump- like distortions. 



features of this type (see, for example, [78|) and hence ob- 
servational constraints for k >, 0.1 Mpc -1 may be weaker 
than what is presented in Fig. [3] Note that in most of 
this nonlinear regime we find essentially no constraint on 
Am-, which is consistent with what would be expected in 
a more comprehensive treatment. 

The IR cascading bump in the primordial power spec- 
trum will be accompanied by a cor resp onding nongaus- 
sian feature in the bispectrum 0, |22| . From the per- 
spective of potentially observing this signal it is most 
interesting if this feature is located on scales probed by 
CMB experiments. (There is also the fascinating possi- 
bility that the nongaussianity from IR cascading could 
show up in the large scale structure as in [7^, H(| El, HI] • 
We leave a detailed discussion to future studies.) To get 
some intuition into what kinds of features in the CMB 
scales are still allowed by the data we focus on an example 
with Am = 2.5 x 10 -10 which, using ([5]), corresponds to a 
reasonable coupling value g 2 ~ 0.01. We take the bump 
to be located at fciR, = 0.01 Mpc -1 and fix the remaining 
model parameters to A s = 2.44 x 10~ 9 , n s — 0.97 (which 
are compatible with the usual values). This sample bump 
in the power spectrum is illustrated in the top panel of 
Fig. [Hand is consistent with the data at 2a. In the bot- 
tom panel of Fig. [4] we plot the associated angular CMB 
TT spectrum. This example represents a surprisingly 
large spectral distortion: the total power in the feature 
as compared to the scale invariant vacuum fluctuations 
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TABLE III: constraints on the standard ("vanilla") cosmolog- 
ical parameters for the multiple-bump model. All error bars 
are 95.4% confidence level. 
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is -Phiimn/fya.n ~ 0.1, evaluated at the peak of the bump. 
In [22| we discuss the nongaussianity associated with this 
feature. 



B. Multiple Bursts of Particle Production 

Next, we consider a slightly more complicated scenario: 
multiple bursts of particle production leading many lo- 
calized features in the power spectrum. For simplicity 
we assume that all bumps have the same magnitude 
At = Am and we further suppose a fixed number of e- 
foldings SN between each burst of particle production. 
This implies that the location of the i-th bump will be 
given by k% = e^~^ A ki where A is a model parameter 
controlling the density of features. We take the num- 
ber of bursts, n, to be sufficiently large that the series 
of features extends over the whole observable range. In 
the next section we will see that these assumptions are 
not restrictive and that many well-motivated models are 
consistent with this simple set-up. 

Our multi-bump model, then, has three parameters: 
Am, k\ and A. We take the prior on the amplitude to 
be Am < 25 x 1CP 10 as in section llV A[ If the features are 
very widely spaced, A > 1, then the constraint on each 
bump will obviously be identical to the results for the 
single-bump case presented in the section IIV Al Hence 
the most interesting case to consider is A <, 1 so that 
the bumps are partially overlapping. Our prior for the 
density of features is therefore < A < 1. Finally, the 
location of the first bump will be a historical accident in 
realistic models, hence we marginalize over all possible 
values of k\ and present our constraints and 2-d likeli- 
hood plots in the space of A\^ and A. This marginal- 
ized likelihood plot is presented in Fig. [5] In table IIIII 
we present the best-fit values for the vanilla cosmological 
parameters. 

From the likelihood plot, Fig. there is evidently a 
preferred value of the feature spacing, roughly A ~ 0.75, 
for which the constraints are weakest. This can be un- 
derstood as follows. For very high density A — > the 
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FIG. 4: The top panel shows a sample bump in the power 
spectrum with amplitude Am, = 2.5 x 1CP 10 which corre- 
sponds to a coupling g 2 ~ 0.01. The feature is located at 
km. ~ 0.01 Mpc" 1 . This example represents a distortion of 
O(10%) as compared to the usual vacuum fluctuations and is 
consistent with the data at 2a. The bottom panel shows the 
CMB angular TT power spectrum for this example, illustrat- 
ing that the distortion shows up mostly in the first peak. 
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FIG. 5: Marginalized posterior likelihood contours for the 
parameters Air and A (the feature amplitude and spacing, 
respectively) of the multiple-bump model. Black and grey re- 
gions correspond to values allowed at 95.4% and 99.7% con- 
fidence levels, respectively. 
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localized features from IR cascading smear together and 
the total power spectrum (g]) is P(k) ~ A s (k/k ) 7ls ^ 1 +C 
where the size of the constant deformation scales linearly 
with the density of features: C oc A -1 . Therefore, the 
upper bound on the amplitude Am, should scale linearly 
with A. Indeed, this linear trend is very evident from 
Fig.[5]in the small-A regime. This linear behaviour must 
break down at some point since as the features become 
infinitely widely spaced the constraint on Am must go 
to zero. This explains the bump in the likelihood plot, 
Fig. El near A - 0.75. 

In passing, notice that the behaviour P(k) ~ 
A s {k/ko) nB ~ 1 +C for A -C 1 also explains why the best- 
fit A s in table IIIII is somewhat lower than the standard 
value and why the spectral tilt n s — 1 is somewhat more 
red. 

To get some intuition for the kinds of multi-bump dis- 
tortions that are allowed by the data, we consider an ex- 
ample with Air = 1 x 10~ 9 , A = 0.75 and fix the vanilla 
parameters to A s — 1.04 x 10~ 9 , n s = 0.93. This choice 
of parameters is consistent with the data at 2a and cor- 
responds to a reasonable coupling g 2 ~ 0.02. In Fig. [S] 
we plot the primordial power spectrum P(k) and also the 
CMB TT angular power spectrum for this example. 
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V. PARTICLE PHYSICS MODELS 

From the low energy perspective one expects interac- 
tions of the type Jl} to be rather generic, hence particle 
production during inflation may be expected in a wide 
variety of models. In this section we consider some ex- 
plicit examples in string theory and SUSY in order to 



FIG. 6: The top panel shows a sample multiple-bump distor- 
tion with amplitude Am. = 1 x 10~ 9 which corresponds to 
g 2 ~ 0.02. The feature spacing is A = 0.75. We take the 
vanilla parameters to be A a = 1.04 x 10 -9 , n s — 0.93 so that 
the scale of inflation is slightly lower than in the standard 
scenario and the spectral tilt is slightly redder. The bottom 
panel shows the CMB angular TT power spectrum for this 
example. 
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show how such models may be obtained microscopically 
and also to provide the proof of concept that realistic 
models do exist where 4>i are in the observable range. 



A. Open String Inflation Models 

String theory inflation models may be divided into two 
classes depending on the origin of the inflaton: closed 
string models and open string models. In the former 
case the inflaton is typically a geometrical modulus as- 
sociated with the compactification manifold (examples 
include racetrack inflation 18311 . Kahler modulus inflation 
[8"j ] and Roulette inflation [85|). In the latter case the 
inflaton is typically the position modulus of some mo- 
bile D-brane 5 living in the compactification manifold (ex- 
amples include brane inflation [13] such as the warped 
KKLMMT model ||, D3/D7 inflation [H] and DBI in- 
flation [90]). In open string inflation models there may 
be, in addition to the mobile inflationary brane, some 
additional "spectator" branes. If the mobile brane col- 
lides with any spectator brane during inflation then some 
of the stretched string states between these branes will 
become massless at the moment when the branes are co- 
incident [3, Ul| > precisely mimicking the interaction {1} . 
Thus, we expect particle production, IR cascading and 
the bump-like features described above to be a reason- 
ably generic prediction of open string inflation. 



B. String Monodromy Models 

A concrete example of the heuristic scenario discussed 
in the last subsection is provided by the brane mon- 
odromy and axion monodromy string theory inflation 
models proposed in [l8l Il9l . [20| . In the original brane 
monodromy model [l8j one considers type IIA string the- 
ory compactified on a nil manifold that is the product of 
two twisted tori. The metric on each of these twisted tori 
has the form 

ds 2 

— = L 2 ui du\ + L 2 U2 du 2 2 + L 2 x {dx + Muxdu 2 ) 2 (6) 

where x' = x — ^U\u 2 and M is an integer flux number. 
The dimcnsionless constants L Ul , L U2 and L x determine 
the size of the compactification. 

Inflation is realized by the motion of a D4-brane along 
the direction u\ of the internal manifold. The D4 spans 
our large 3-dimensions and wraps a I-cycle along the di- 
rection u 2 of the internal space. The size of this 1-cycle, 
in string units, is given by 



T 2 



L 2 x M 2 u\ 



(7) 



5 One notable exception is inflation driven by the open string 
tachyon, for example nonlocal string field theory models |8b|. 



Hence, the brane prefers to minimize its world-volume 
by moving to the location u\ = where this I-cycle has 
minimal size. This preference gives a potential to the 
D4-brane position which goes like V ~ u\ in the large ui 
regime that is relevant for large field inflation. 

In [TJ| it was shown that this scenario allows for the 
inclusion of a number of spectator branes stabilized at po- 
sitions Mi = i/M (with i integer) along the inflationary 
trajectory. As the mobile inflationary D4 rolls through 
these points particles (strings) will be produced and the 
resulting distribution of features will look nearly identi- 
cal to the simple multi-bump scenario studied in section 
IIVBI To see this, let us now determine the distribution 
of features that is predicted from brane monodromy in- 
flation. The canonical inflaton can be related to the 
position of the mobile D4 as 



= Bu{ 



(8) 



where B is a constant with dimensions of mass that de- 
pends on model parameters. Hence, the effective poten- 
tial during inflation has the power-law form 



V(<P) = A P P 



(9) 



For the simplest scenario described above one has p — 
2/3. However, the formulas (|8l9p still hold for the variant 
considered in [1 81 ] with p — 2/5 as long as one replaces u\ 
by a more complicated linear combination of coordinates. 
These relations also hold for axion monodromy models 
[lj| with p — 1 and u\ replaced by the axion, c, arising 
from a 2-form RR potential integrated over a 2- 

cycle E 2 - For all models of the form the number 
of e-foldings N from <f> = <f>(N) to the end of inflation 
(which occurs at <j) = pj \[2 when the slow roll parameter 
e(</>) = 1) is given by 



N = 



1 4> 2 (N) 



M 2 



2p 
2PM 2 '' 



P 
4 



,2/p 



P 

4 



(10) 



Since the spectator branes are located at u\ — i/M the 
bursts of particle production must occur at times N = Ni 
where 



Ni 



fyM 2 



2/p 



p 

4 



(11) 



The location k = ki of the i-th feature is defined, roughly, 
by the scale leaving the horizon at the moment N — 
Ni. Hence, the distribution of features for brane/axion 
monodromy models is given by 



hi 



S B 2 i 2 /p - 1 
4 



(12) 



with p = 2/3 or p = 2/5 for brane monodromy and p = 1 
for axion monodromy. In (|12[) the dimensionless number 
B depends on model parameters. 
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Although the distribution of features (|12[) is not ex- 
actly the same as the evenly space distribution considered 
subsection liVB| the two are essentially indistinguishable 
over the range of scales which are probed observation- 
ally (corresponding to roughly 10 e- foldings of inflation). 
The reason for this is simple: the inflaton is nearly con- 
stant during the first 10 e-foldings of inflation and hence 
8N ~ 5<j) ~ 8u\ within the observable region. It follows 
that ki = e^ 1 ' k% to very good approximation for a 
huge class of models. To see this more concretely in the 
case at hand, let us compute dN/dui from (fT0"|) . It is 
straightforward to show that 



dN 1 1 



du i P p [2e(4>)\ 



!-p/2 iM 



B V 



where 



(13) 



(14) 



is the usual slow roll parameter. Observational con- 
straints on the running of the spectral index imply 
that e(</>) cannot change much over the observable 10 e- 
foldings of inflation. Since dN/dui = const to very high 
accuracy it follows trivially that N = N(u\) is very close 
to linear and ki = e' t_1 ^ A fci as desired. 

In the context of axion monodromy inflation models 
[l^ | the multiple bump features discussed here will be 
complimentary to the oscillatory features described in 
(20| which result from the sinusoidal modulation of the 
inflaton potential by instanton effects. If the bursts of 
particle production are sufficiently densely spaced, then 
signal from IR cascading may appear oscillatory, how- 
ever, it differs from the effect discussed in [2(| in both 
physical origin and also in functional form. 

Let us now estimate the effective value of the couplings 
g\ appearing in the prototype interaction ([3]) that are 
predicted from the simplest brane monodromy model. 
A complete calculation would involve dimensionally re- 
ducing the DBI action describing the brane motion and 
requires knowledge of the full 10-dimensional geometry 
with the various embedded branes. For our purposes, 
however, a simple heuristic estimate for the collision of 
two D4-branes will suffice. When N D-branes become 
coincident the symmetry is enhanced from U(1) N to a 
U(N) Yang Mill gauge theory. The gauge coupling for 
this Yang Mills theory is given by 



2 

9ym 



5s (2tt) 



(15) 



where L is the volume of the 1-cycle that the D4 branes 
wrap and is given by ([?]). If the inflationary brane is at 
position u% and the i-th spectator brane is at u\ t i then 
the distance between the two branes is given by 



d 2 = a 1 L 2 ai (ux - ui^Y 



(16) 



The mass of the gauge bosons corresponding to the en- 
hanced symmetry is 



Mf 



2 

3ym 



A2 



L 2 2 



LlAPu 2 



(17) 



To put this in the prototype form M 2 — g 2 (<fi — (f>i) 2 we 



must first convert to the canonical variable 
formula © with p = 2/3 and 



D 



M 1 ' 2 L ni Lx 



1/2 



67r 2 \Jg s a' 



using the 



(18) 



Next, we must Taylor expand the resulting equation 
about the minimum d> — 6, ■ . We find 



Mi S 



16<? 2 7r 4 1 



ML x Ul,i ^ L 2 2 +L2 M 2 U 2. 

16^ 4 1 



L. r i 



(19) 



(20) 



where on the second line of (f2"0"]) we have used the fact 
that iti.i = i/M (with i integer) in the simplest models. 
We see that the effective couplings gf become larger as 
the D4 unwinds during inflation. (The apparent diver- 
gence for ui : i — in the formula (|20[) is an artifact of 
the fact that the relation ([8]) is not valid at small values 
of Mi. This will not concern us here since inflation has 
already terminated at the point that our formulas break 
down.) 

To compute the amplitude of the bump-like feature 
produced by brane monodromy inflation we should take 
into account also combinatorial factors. When two 
branes become coincident the symmetry is enhanced from 
U(l) 2 to U(2) so there are 2 2 — 2 = 2 additional mass- 
less spin-1 fields appearing at the brane collision. Thus, 
using equation ([5]), the amplitude of the feature that will 
be imprinted in the CMB is 



Ai, eS = 2 x (2 2 - 2) x 



1.01 • 10 



-e 1B /4 



(21) 



where the extra factor of 2 counts the polarizations of the 
massless spin-1 fields. This combinatorial enhancement 
can be much larger if the inflationary brane collides with 
a stack of spectators. 

The above discussion is predicated on the assumption 
that the original brane monodromy set-up [18j is supple- 
mented by additional spectator branes. This may seem 
like an unnecessary contrivance, however, in order for this 
model to reheat successfully it may be necessary to in- 
clude spectator branes. For example, with the reheating 
mechanism proposed in 91( semi-realistic particle phe- 
nomenology can be obtained by confining the standard 
model (SM) to a D6 brane which wraps the compact 
space. In order to reheat into SM degrees of freedom we 
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orient this brane so that its world-volume is parallel to 
the mobile (inflationary) D4. In this case the end of in- 
flation involves multiple oscillations of the D4 about the 
minimum of its potential. At each oscillation the D4 col- 
lides with the D6 and SM particles are produced by para- 
metric resonance preheating [HUG^l- However, due to the 
periodic structure of the compactification, D4/D6 colli- 
sions will necessarily occur also during inflation, leading 
to IR cascading features in the CMB. 

The timing of these D4/D6 collisions was computed in 
[9l| for the minimal p — 2/3 brane monodromy model, 
assuming the same choices of parameters used in [18| . 
For this particular case there is only one collision (and 
hence one feature) during the first 10 e-foldings of infla- 
tion and the phenomenology is essentially the same as 
that considered in subsection IIV Al What is the ampli- 
tude of this feature? Assuming, again, the parameters 
employed in fl^| and noting that the first collision takes 
place at i — 13 [9l| equation (f2T)]) gives g\ = 0.001. From 
(l2Tj) we find the effective amplitude of the feature to be 
Ai jB ff/A s = 0.01. This value is well within the observa- 
tional bounds derived in subsection IIV Al 

We stress that the conclusions in the previous para- 
graph apply only for the particular choice of model pa- 
rameters employed in [l8|. There exist other consistent 
parameter choices for which the simplest brane mon- 
odromy model predicts a much higher density of features 
with much larger amplitude. 

Note that both brane and axion monodromy models 
may be used to realize trapped inflation [15| . Here we 
are restricting ourselves to the large-field regime where 
the potential V = y^~ p (\P is flat enough to drive infla- 
tion without the need for trapping effects. For a given 
choice of parameters one should verify that this classical 
potential dominates over the quantum corrections from 
particle production. 



C. A Supersymmetric Model 

Another microscopic realization of multiple bursts of 
particle production and IR cascading during inflation 
which does not rely on string theory can be obtained 
from the so-called "distributed mass" model derived in 
[2ll ] with warm inflation [l7j in mind, however, the theory 
works equally well for our scenario. This model is based 
on M = 1 global SUSY and allows for the inclusion of 
multiple points along the inflationary trajectory where 
both scalar degrees of freedom and also their associated 
fermion superparteners become massless. The distribu- 
tion of features in this set-up is essentially arbitrary. 



VI. CONCLUSIONS 

In this paper we have studied the observational con- 
straints on models with particle production during infla- 
tion. We have focused on the simple prototype model ([T]) 



for each burst of particle production, however, we expect 
that our qualitative results will apply also to more com- 
plicated models (for example with gauged interactions or 
fermion iso-inflaton fields) and perhaps also to the case 
of inflationary phase transitions. We find no evidence for 
a detection of the features associated with particle pro- 
duction and IR cascading, however, it is interesting to 
note that rather large localized features are still compat- 
ible with the data. Our results differ significantly from 
previous studies as a result of a more realistic treatment 
of the cosmological perturbations in models with parti- 
cle production. The bounds we have derived on g 2 will 
play a crucial role in assessing the detectability of the 
nongaussianity produced by particle production and IR 
cascading. 

We have also discussed the implications of our results 
for popular brane/axion monodromy string theory infla- 
tion models. Successful reheating in these constructions 
may require the inclusion of spectator branes which col- 
lide with the mobile D4-brane during inflation and hence 
we expect CMB features to be a fairly generic predic- 
tion. We have shown that brane/axion monodromy mod- 
els predict a distribution of bump-like features which are 
evenly spaced in lnfc over the observable range of scales. 
In the case of axion monodromy this multiple-bump spec- 
tral distortion is complimentary to the oscillatory fea- 
tures discussed in [2(| ■ We have also estimated the mag- 
nitude of these bump-like features in terms of model pa- 
rameters. 

One motivation for the present study was to determine 
the extent to which microscopically realistic models such 
as (Tl]) can reproduce the localized "glitches" in the power 
spectrum that have been detected (albeit with marginal 
significance) by several previous studies. These anoma- 
lies can be classified as follows: 

1. Localized power excesses: 

In both and [33| power spectra with local- 
ized spikes were studied and in both cases marginal 
evidence was found for a detection of such fea- 
tures. In [3~H a non-parametric reconstruction of 
the power spectrum was performed and the result 
is marginally consistent with a power law every- 
where, however, several localized spikes are evident 
in the reconstruction. 

Localized excesses are naturally obtained in our 
model (fT]). Sadly, however, we did not find that 
our model fits the data significantly better than 
the simplest slow roll inflation scenario. This does 
not necessarily imply a disagreement with 0, HH 
since we use a different shaped feature and different 
data sets. (Indeed, when the authors of [33[ repeat 
their analysis using the WMAP 5-year data they 
do not obtain a detection [92| . consistent with our 
findings.) 

2. Localized power deficits: 
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In [24[ the Richardson-Lucy deconvolution algo- 
rithm was used to perform a non-parametric recon- 
struction of the primordial power spectrum which 
displayed a prominent IR cut-off near the horizon. 
In [3_4j a similar analysis was performed and the 
reconstructed power spectrum displays a localized 
dip in power near k ~ 0.002 Mpc -1 . 

Localized deficits can be produced by our model 
([3]) but only in a rather contrived way. Hence, we 
have not focused on such features in section [TV] 

3. Damped oscillations: 

In [ll|, [l^, HE HE] power spectra with superim- 
posed ringing patterns were studied. Such features 
provide a marginally improved fit over the simplest 
power-law model. 

As we have discussed in the introduction, damped 
oscillatory "ringing" features are not predicted by 
inflationary particle production. Nor is it clear if 
such features are predicted by models with phase 
transitions. (Of course damped oscillations can be 
obtained from a toy model with a step in V(<fi). 
However, it may be difficult to obtain such a po- 



tential from realistic micro-physics; generically one 
expects that any sharp features in V(</>) will be 
smoothed out by quantum corrections.) 

Finally, let us note that features of the type studied 
here will lead to other observables beyond the distortion 
of the primordial power spectrum. In particular, bumps 
in P{k) will lead to features in the tensor spectrum (re- 
sulting from the sourcing of gravitational waves by scalar 
fluctuations at second order in perturbation theory) and 
also, possibly, black hole production. In [93l |94| these 
effects were estimated assuming a power spectrum which 
is qualitatively similar to ours. As discussed in [5j, infla- 
tionary particle production will also lead to potentially 
large localized nongaussian features in the bispectrum 
(and higher order statistics) of the cosmological fluctua- 
tions. These nongaussianities will be discussed in detail 
in an upcoming work [22l |. 
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